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Abstract 

Determination of the neutrino mass hierarchy in a reactor neutrino experiment 
at medium baseline is discussed. Observation of the interference effects between 
the Am^ and Am| 2 oscillations enables a relative measurement independent of 
the knowledge of the absolute mass-squared difference. With a 20 kton liquid scin- 
tillator detector of 3%/y/E(MeV) energy resolution, the Daya Bay II Experiment 
at a baseline of ~ 50 km from reactors of total thermal power 36 GW can determine 
the mass hierarchy at a confidence level of Ax^ H ~ (10 -j- 12) (3 -j- 3.5c) in 6 years 
after taking into account the real spatial distribution of reactor cores. We show 
that the unknown residual energy non-linearity of the liquid scintillator detector 
has limited impact on the sensitivity due to the self-calibration of small oscillation 
peaks. Furthermore, an extra increase of A%mh — ^ (9) can De obtained, by in- 
cluding the precise measurement of the effective mass-squared difference Am^ of 
expected relative error 1.5% (1%) from ongoing long-baseline muon neutrino dis- 
appearance experiments. The sensitivities from the interference and from absolute 
measurements can be cross checked. When combining these two, the mass hierarchy 
can be determined at a confidence level of Axjyj H ~ (15 -j- 20) (4(7 ) in 6 years. 
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1 Introduction 



An unexpectedly large value of neutrino mixing angle #13 was measured by the Daya Bay 
Reactor Neutrino Experiment (DYB) [HI2], together with other consistent evidences from 
the reactor [3J[I] and accelerator [SJE] neutrino experiments. It opens a gateway to the 
measurement of neutrino mass hierarchy (i.e., the sign of Am^ or Am§ 2 ) and the leptonic 
CP- violating phase (5cp). A large value of #13 makes both above measurements easier 
in the next generation of neutrino oscillation experiments. Possible information on the 
neutrino mass hierarchy (MH) can be extracted not only from the matter-induced oscil- 
lations of long-baseline accelerator neutrino experiments [5HTUj and atmospheric neutrino 
experiments [TT|[T2]. but also from the vacuum oscillation in reactor neutrino experiments 
at medium baseline [TBTI24"] . 

The MH sensitivity of a future reactor neutrino experiment comes from the interfer- 
ence effect of two separated oscillation modes [H)J[T7] (i.e., two fast oscillations driven by 
Am^ and Am^). The relative sizes of [Am^l and |Am 2 2 | and the non- maximal value 
of 812 make it possible to determine the MH in vacuum, immune from the uncertainty 
of the Earth density profile and ambiguity of the CP- violating phase in atmospheric and 
long-baseline accelerator neutrino oscillation experiments. Such a measurement is chal- 
lenging. The energy resolution as the size of Am^/I Am^l and event number of several 
tens of thousands are the minimal requirements [17]. Moreover, the spatial distribution 
of reactor cores [20H22] and non-linearity of the energy response [21] may also degrade 
the sensitivity. 

Besides the interference effect in vacuum oscillations, precision measurement of the 
flavor-dependent mass-squared differences (i.e., Am^ e , Am^ and Am^ T ) [25l|26] can also 
contribute to the MH sensitivity. Am^ e from reactor experiments and Am^ from long- 
baseline accelerator muon-neutrino disappearance experiments are different combinations 
of Amj! and Am| 2 - Therefore, a comparison of the effective mass-squared differences with 
different flavors can discriminate the neutrino MH if they are measured with sufficient 
precision. Considering the MH determination from the Daya Bay II reactor neutrino 
Experiment, the sensitivity from the interference effect can be improved by including 
an accurate measurement of the effective mass-squared difference Am^ from accelerator 
neutrino experiments. 

We first give a brief description on the effective mass-squared differences in different 
oscillation channels in Section 2. Statistical analysis of Daya Bay II is introduced in 
Section 3. Impact from the energy non-linearity of the liquid scintillator (LS) detector is 
discussed in Section 4. Improved MH sensitivity with the absolute mass-squared difference 
Am 2 e and external Am^ measurements is presented in Section 5. Finally, we conclude 
in Section 6. 



2 Effective Mass- Squared Differences 

In the standard three-neutrino mixing scheme the survival probability for the a-flavor 
neutrinos is given by [7] 

P{y a v a ) = P(P a V a ) = 1 - 4|[/ Q3 | 2 |f/ Ql | 2 sin 2 A 31 

- 4|^3| 2 |^2| 2 sin 2 A 32 (1) 
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- 4|f/ a2 | 2 |[4 1 | 2 sin 2 A 2 i, 

where A^ = AmfjL/AE, Am 2 ,- = mf — mh and U a i is the element of the leptonic mixing 
matrix. Note that only two of the three A^ are independent because we have the relation 
Am§! = Am| 2 + Am^ from their definitions. 

Strong hierarchy between the magnitudes of Am^ and Am^ (or Am§ 2 ) is achieved 
from the analysis of solar (or KamLAND) and atmospheric (or long-baseline accelerator) 
neutrino oscillation data [2TH29] . At baselines that oscillation driven by the atmospheric 
mass-squared difference dominates, we have A 2 i <C 1 and the corresponding oscillation 
terms are extremely small. In such a neutrino oscillation experiment, it is impossible to 
distinguish between the two neutrino mass hierarchies because two degenerate solutions 
with different hierarchies can generate the identical neutrino energy spectrum. The ab- 
solute values of Am^ (or Am 2 2 ) in the two solutions are different due to the non-zero 
and positive value of Am^. However, we can define an effective mass-squared difference 
as the linear combination of Am^ and Am§ 2 : 

Am 2 Q = ijaAmj! + (1 - ?7 a )Am 2 2 

= Am 2 2 + 7] a Am 2 2l = Am 2 x - (1 - 7] a )Am 2 2l , (2) 

where the absolute value of the effective mass-squared difference keeps invariant for both 
degenerate solutions. Generally, the value of r] a varies for different oscillation channels 
because of the flavor-dependent amplitudes in these oscillation probabilities. Using the 
power series of A 21 , we can rewrite the probability in Eq. ([T|) as 

P{v a -> v a ) = 1 - 4|[/ a3 | 2 (l - |£4 3 | 2 ) sin 2 A aa 

- ±\U a3 \ 2 [\U al \ 2 (l-r) a )-\U a2 \ 2 r) a ] A 21 sin2A aa (3) 

- 4|f/ Q2 | 2 |[/ al | 2 sin 2 A 21 -(9(A 2 1 ), 

where A aa = Am^ a L/4E. Therefore, we obtain an effective two-neutrino oscillation 
scheme up to the linear term of A 2 i, if the parameter rj a satisfy the following relation 





\U a l\ 


to 


|£4i| 


2 + 


\Ua2 


2 



(4) 



Using the standard parametrization of the leptonic mixing matrix [7], we get the effective 
mass-squared differences in Eq. (j2J) for different channels of neutrino oscillations 

Am 2 e ~ cos 2 6'i 2 Am 2 1 + sin 2 8 12 Am 2 2 , (5) 
Am^ ~ sin 2 9± 2 Am 31 + cos 2 9 12 Am 2 2 + sin 29± 2 sin 6*13 tan 9 23 cos 5Am 21 , (6) 
Am 2 T ~ sin 2 12 Am 2 1 + cos 2 6 l i 2 Am 2 2 — sin26'i 2 sin 9i 3 cot 6* 2 3 cos <5Am 21 , (7) 

where terms at the order of (9(sin 2 9i 3 Am\ 1 ) have been neglected for simplicity. We can 
also calculate the differences of the effective quantities between different flavors as 

|Am 2 e | — |Am 2 M | = ±Am 21 (cos 2#i 2 — sin 29\ 2 sin #13 tan 6* 2 3 cos 5) , (8) 
I Am 2 J — |Am 2 r | = ±2Am 21 sin 2^ 12 sin 13 esc 2^23 cos 5 , (9) 
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where the positive and negative signs correspond to normal and inverted mass hierar- 
chies, respectively. Therefore, comparisons of oscillations between different flavors may 
distinguish the MH and even tell us possible information on the CP-violating phase. 

One must keep in mind that the effective mass-squared differences defined here are 
only MH-invariant under the condition of A 2 i <C 1. In the reactor neutrino experiment 
at a medium baseline (A21 ~ 1), where all the oscillation modes and their interference 
terms are measurable, the absolute value of Am^ is not invariant by changing the sign 
of the neutrino MH. However, it is still close to Am^, rather than Am^ and we can get 
an additional MH sensitivity with the inclusion of a prior Am^ measurement. 

For a reactor neutrino experiment at medium baseline, corrections to the mass-squared 
differences from the terrestrial matter effect are around 1% and the induced uncertainties 
are negligibly small (less than 0.1%). On the other hand, in the muon- neutrino disappear- 
ance channel of long-baseline accelerator neutrino experiments, the matter corrections are 
suppressed by the smallness of #13 and only at the level of 0.2% for the baselines of several 
hundreds kilometers (e.g., 295 km for T2K [3U] and 735 km for NOvA [3T]). Moreover, 
the different signs in the matter potentials of neutrino and antineutrino oscillations are 
also favorable to increase the discrepancy of different mass-squared differences. 



3 Statistical Analysis 

The 20 kt liquid scintillator (LS) detector of Daya Bay II Experiment [20TT22] will be 
located at equal baselines of 52 km away from two reactor complexes (36 GW in total). 
In this study we use nominal running time of six years, 300 effective days per year, and a 
detector energy resolution 3%/ v /£(MeV) ClS db benchmark. A normal MH is assumed to 
be the true one (otherwise mentioned explicitly) while the conclusion won't be changed 
for the other assumption. The relevant oscillation parameters are taken from the latest 
global analysis [27] as Am 2 21 = 7.54 x lO^eV -2 , (Am^ + Am| 2 )/2 = 2.43 x lO^eV -2 , 
sin 2 6*i3 = 0.024 and sin 2 6*12 = 0.307. The CP-violating phase will be specified when 
needed. Finally, the reactor antineutrino flux model from Vogel et al. [32] is adopted in 
our simulation.]. 

To obtain the sensitivity of the proposed experiment, we employ the least squares 
method and construct a standard x 2 function as following: 

^ [M t -T t (l + E fc ^ fc e fc )] 2 4 

i=l k K 

where M, is the measured neutrino events in the z-th energy bin, Tj is the predicted 
reactor antineutrino flux with oscillations, is the systematic uncertainty, is the 
corresponding pull parameter, and is the fraction of neutrino event contribution of 
the fc-th pull parameter to the i-th energy bin. The considered systematic uncertainties 
include the correlated (absolute) reactor uncertainty (2%), the uncorrelated (relative) 
reactor uncertainty (0.8%), the flux spectrum uncertainty (1%) and the detector-related 

1 We have tried both the calculated [32] and the new evaluations [331133] of the reactor antineutrino 
fluxes. The discrepancy only influences the measurement of 612. Both evaluations give consistent results 
on the MH determination. 
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Figure 1: The MH discrimination ability for the proposed reactor neutrino experiment 
as functions of the baseline (left panel) and the detector energy resolution (right panel) 
with the method of the least squares function in Eq. (I10|) . 
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Table 1: Summary of the power and baseline distribution for the Yangjiang (YJ) and 
Taishan (TS) reactor complexes, as well as the remote reactors of Daya Bay (DYB) and 
Huizhou (HZ). 

uncertainty (1%). We use 200 equal-size bins for the incoming neutrino energy between 
1.8 MeV and 8.0 MeV. 

We can fit both the normal MH and inverted MH with the least squares method 
and take the difference of the minima as a measurement of the MH sensitivity. The 
discriminator of the neutrino MH can be defined as 

AX^H=lxLn(N)-xLn(I)l, (11) 

where the minimization process is implemented for all the relevant oscillation parameters. 
Note that two local minima for each MH [Xmin(N) and Xmm(-0] can be located at different 
positions of |A?rig e |. This particular discriminator is used to obtain the optimal baseline 
and to explore the impact of the energy resolution, which are shown in the left and right 
panels of Figure 1. Ideally a sensitivity of ^Xmu — 16 can be obtained at the baseline 
around 50 km and with a detector energy resolution of 3%. 

The baselines to two reactor complexes should be equal. The impact of unequal 
baselines is shown in the left panel of Figure 2, by keeping the baseline of one reactor 
unchanged and varying that of another. A rapid oscillatory behavior is observed and 
demonstrates the importance of baseline differences for the reactor cores. To evaluate 
the impact from the spacial distribution of individual cores, we take the actual power 
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Figure 2: The variation (left panel) of the MH sensitivity as a function of the baseline 
difference of two reactors and the comparison (right panel) of the MH sensitivity for the 
ideal and actual distributions of the reactor cores. 
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Figure 3: Two classes of typical examples for the residual non-linear functions in our 
simulation. 

and baseline distribution of each core of the Yangjiang (YJ) and Taishan (TS) nuclear 
power plant, shown in Table 1. The remote reactors in the Daya Bay (DYB) and the 
possible Huizhou (HZ) power plant are also included. The reduction of sensitivity due to 
the actual distribution of reactor cores is shown in the right panel of Figure 2, which gives 
a degradation of ^Xmh — 5- in all the following studies, the actual spacial distribution 
of reactor cores for the Daya Bay II Experiment is taken into account. 



4 Energy Non-Linearity Effect 

The detector energy response is also crucial for Daya Bay II since a precise energy spec- 
trum of reactor neutrinos is required. Assuming the energy non-linearity correction is 
imperfect, we study its impact to the sensitivity by including in our simulation a residual 
non-linearity between the measured and expected neutrino spectra. Assume the detector 
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energy non-linearity has the form as 



-Etruc 1 + Pi exp(-p 2 E tIUC ) ' 



where E rec and E true are the reconstructed and true kinematic energy of the positron from 
the inverse beta decay, respectively. The parameter of p , pi and p 2 describe the shape and 
magnitude of the non-linear functions. We assume that, after the non-linearity correction, 
the residual non-linearity in measured energy spectrum also has the same function form. 
The conclusion won't be changed for other residual non-linearity assumptions, because 
we will use a quadratic function in the predicted spectrum, different from the measured 
spectrum. We fix p 2 = 0.2 /MeV and vary po and p\ as 

Po = sign x sizei and pi = sign x size2 , (13) 

where sign = ±1 determines the slope, sizei and size2 can be a few percent to indicate 
the magnitudes of the residual non-linearity. Two typical examples with sign = ±1 are 
shown in Figure 3. 

By including the residual non-linearity in M; of Eq. (fTUj) , we obtain in Figure 4 updates 
on the distributions of the A% 2 function for both true and false neutrino MH, where 
normal (inverted) MH is assumed to be the true one in upper (lower) panels. Different 
classes of the non-linear functions may induce different effects on the MH determination. 
Comparing to that without non-linearity effects as shown in the right panel of Figure 2, 
the non-linearity with positive sign (sign = +1) will increase the discrepancy between 
two neutrino MH scenarios for the normal true MH, but decrease the discrepancy for the 
inverted true MH. The non-linear functions with negative sign (sign = — 1) have opposite 
effect. Only when the size of non-linearity is as small as 0.5% this effect can be ignored. 

For the reactor antineutrino experiment at intermediate baselines (~ 50 km), we can 
observe multiple peaks of the Am 2 ^ induced oscillation. Each peak position carries the 
information of the same Am 2 e . This redundancy can be used to evaluate the energy scale 
at different energies. Therefore, we can measure to some extent the energy non-linearity 
by the spectrum itself. We call this effect as the self-calibration of the spectrum. To 
illustrate, we consider a test quadratic non-linear function in the fitting process, 

-p^ ~ 1 + g + <?i£truc + q 2 El ue , (14) 

-'-'true 

where the central values of three parameters are arbitrary, but the uncertainties are as- 
sumed to be the same size as the function in Eq. ( TT2"j) . We can define the new least-square 
function by using the test non-linearity function in Tj and including the corresponding 
pull terms Xnl = ^A^*?*) 2 anc ^ derive the MH sensitivity taking into account the 
non-linearity and self-calibration effects. Considering different sign and size of nonlinear- 
ity defined in Eq. (13), we illustrate the A% 2 function with the self-calibration in fitting 
in Figure 5, where the normal (inverted) MH is assumed to be the true one in the upper 
(lower) panels. First we can observe that the degeneracy ambiguity induced by the non- 
linearity effect can be removed by fitting the parameters of the test non-linearity function 
and both classes of non-linear functions give the consistent sensitivity of MH determina- 
tion (see the left and right panels of Figure 4 and 5). Tiny differences can be noticed for 
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Figure 4: Effects of two classes of energy non-linearity models in the determination of 
the neutrino MH without the self-calibration in fitting. The normal (inverted) MH is 
assumed to be the true one in the upper (lower) panels. The sign and size of the non- 
linear parameters are indicated in the legend. 
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Figure 5: Effects of two classes of energy non-linearity models in the determination of the 
neutrino MH with the self-calibration in fitting. The normal (inverted) MH is assumed 
to be the true one in the upper (lower) panels. The sign and size of the non-linear 
parameters are indicated in the legend. 
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different size of the non-linearity because the test quadratic function cannot describe the 
true residual non-linearity accurately. Second, the width of the A% 2 functions in Figure 
5 is broadened compared to Figure 4. This is because additional uncertainties from the 
non-linearity parameters are introduced, which can be translated to the uncertainty of 
the neutrino spectrum and finally to the accuracy of the oscillation parameters. 



5 Improvement with external ArajL measurements 

Due to the difference between Am^ and Am^, precise measurements of the mass-squared 
differences can provide additional sensitivity to MH, besides the sensitivity from the 
interference effects. 

To incorporate the contribution from the Am^ measurement in the long-baseline 
muon-neutrino oscillation experiments, we define the following the extra pull function 



(|Am 2 I - |Am 2 ,|) 2 
XUI*<I) = ' <15) 



where Am^ and a(Am^) are the central value and la uncertainty of the measurement 
respectively. The combined x 2 function is defined as 

X 2 all = Xrea + Xpuil ( I Am L I ) • ( 1 6 ) 

Because two of the three mass-squared differences (Am^, Am^ and Am^) are inde- 
pendent, we choose Am^ and Am^ defined in Eq. ([5]) as the free parameters in this 
work. Proper values of Am^ can be calculated by the relations in Eqs. (8) and (9). 
In practice we need to consider the uncertainties of other oscillation parameters, but the 
CP-violating phase S is almost non-constrained and we can absorb the uncertainties of 
other parameters in that of the phase 6 and consider the whole range of CP-violating 
phase from 0° to 360°. 

As shown in Eq. ([8]), because the relative size of Am^ and Am^ is different for the 
normal and inverted neutrino MH, the extra pull function Xp U \i * n Eq- (fl5l) can give a 
non-zero contribution to the discriminator Ax^ H at the magnitude of 

2 (M rr\ [2 x Am^ (cos 2fli2 - sin2fl 12 sinfl 13 tan fl 23 cos £)] 2 
X P uiil M±1 J~ ^ 2 (Am^) ' [U) 

and correspondingly improve the MH sensitivity for the reactor neutrino experiment at 
medium baselines. 

To illustrate the effect of the external Am^ measurement, we first fix the non-linearity 
with sign = +1, sizei = 2% and size2 = 4%, choose 5 to be 90°/270° (cos 5 = 0) and give 
the separated and combined distributions of the \ 2 functions in Eqs. ([TO]) and ffT6]) in 
Figure 6, where a 1% (left panel) or 1.5% (right panel) relative error of Am^ is assumed. 
The black and red lines are for the true (normal) and false (inverted) neutrino MH, re- 
spectively. The dashed and solid lines are for the reactor-only [in Eq. (fTOj) ] and combined 
distributions. We can get a value of Ax^ H ~ (10 -4- 11) for the reactor-only analysis in 
the least-squares method. As for the contribution from the external Am^ measurement, 
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Figure 6: the reactor-only (dashed) and combined (solid) distributions of the A% 2 function 
in Eq. fllOp and Eq. (j!6p . where a 1% (left panel) or 1.5% (right panel) relative error of 
AwL is assumed and the CP- violating phase (5) is assigned to be 90°/270° (cos 5 = 0) 
for illustration. The black and red lines are for the true (normal) and false (inverted) 
neutrino MH, respectively. The non-linearity in Eq. (112]) is assigned with sign = +1, 
sizei = 2% and size 2 = 4%. 



it is almost negligible if we choose the true MH in the fitting program. However, if the 
fitting MH is different from the true one, the central value of Am^ in the Xp U n function 
will change by two units of the difference in Eq. (]8]), which accordingly results in a sig- 
nificant contribution to the combined \ 2 function. Finally we can achieve A-Xym — 19 
and Ax^ H ~ 14 for the 1% and 1.5% relative errors of the Am^ measurement. 

Next we can discuss the ambiguity of the unknown CP-violating phase 5 and evolution 
of the MH sensitivity with respect to changes of the AmJ error. The Ax^ H dependence 
on different input errors is shown in Figure 7, where the blue, black and red lines stands 
for different values of the CP-violating phase (5 = 0°, 5 = 90°/270° and 5 = 180° 
respectively). In Figure 7, we can notice that the improvement are obvious for an external 
Am^ measurement better than 2% and becomes significant if we can get to the 1% 
level. For the effect of the CP-violating phase, it is most favorable for the value close 
to 180°. The cases of maximal CP violation are in the middle region which are just 
the cases discussed in Figure 6. The ambiguity of the CP- violating phase can induce 
an uncertainty of A^mh - 2 ( 4 ) at ^(AmJJ/IAmJJ ~ 1.5% (1%). The effect of the 
external Am^ measurement can also be viewed as a probe of the CP-violating phase. If 
the improvement is much better than the discussion in Figure 6, a preference of 5 close 
to 180° can be achieved. Otherwise, we may get a nearly vanishing CP- violating phase if 
the situation is totally opposite. 

Current best measurement for Am^ from the MINOS experiment [SB] gives an error 
of 4%. Two new experiments T2K [3U] and NOvA [31] are in operation or construction 
and each of them can reach 1.5% by 2020 after finishing of their nominal running plans 
(5 years of v mode at 750 kW for T2K and 3 years of v mode plus 3 years of v mode 
at 700 kW for NOvA). If these experiments could extend to another 5- year running, it 
might be possible to obtain the precision of 1% which will be useful for the measurement 
of the precision reactor neutrino experiment. 
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Figure 7: The ^Xmk dependence on different input errors of Am^ is illustrated. The 
blue dashed, black solid and red dotted lines stands for different CP values (5 = 0°, 
5 = 90°/270° and 5 = 180° respectively). 



6 Conclusion 

In this work, we have discussed the determination of the neutrino MH using the Daya Bay 
II reactor neutrino experiment at a medium baseline around 50 km away from reactors. 
Precision measurements of the reactor antineutrino spectrum can probe the interference 
effect of two fast oscillation modes (i.e., oscillations induced by Am^ and Am^) and 
sensitive to the neutrino MH. The corresponding sensitivity depends strongly on the 
size of #13, the energy resolution, the baseline differences and energy response functions. 
Moreover, the MH sensitivity can be improved by including a measurement of the effective 
mass-squared difference in the long-baseline muon-neutrino disappearance experiment 
due to flavor dependence of the effective mass-squared differences. 

We have calculated the MH sensitivity taking into account the real spatial distribution 
of reactor complexes, and demonstrated that the residual energy non-linearity of the liquid 
scintillator detector has limited impacts on the sensitivity due to the self-calibration of 
small oscillation peaks. We numerically calculated the sensitivity by assuming two typical 
classes of energy non-linearity functions (2%) and discussed the improvement with the 
external Am^ measurement quantitatively. To conclude, the Daya Bay II Experiment 
could determine the mass hierarchy unambiguously with a confidence level of Axlm ~ 14 
(3.7 a) or Ax^h ~ 19 (4.4 a) in 6 years, for Am^ uncertainty of 1.5% or 1%. 
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